202309.00226v1 


chinaXiv 


第 46 卷 第 9 期 
2023 年 9 月 


T2 MUR 


ARID LAND GEOGRAPHY 


ChinaXiv 合 作 期 刊 


Vol.46 No.9 
Sep. 2023 


不 同 沉陷 应 力 区 土壤 水 分 和 溶质 运 移 的 模拟 试验 


温 K”, 


S er 2 
沿海 丽 i 


车 建 伟 *， ARA. MET 


C1. 中 国 矿 业 大 学 (北京 ) , 煤 让 资源 与 安全 开采 国家 重点 实验 室 ,北京 ”100083; 


2. 内 蒙古 科技 大 学 矿业 与 煤炭 学 院 , 内 蒙古 包头 ”014010) 


摘 要 : 中 国 山西 - 陕 


-内 蒙古 交界 地 开采 沉陷 对 地 表土 壤 扰 动 强 烈 ,引起 地 表 生 态 损 伤 。 为 了 


揭示 开采 沉陷 对 土壤 水 盐 运 移 的 影响 机 制 ,通过 建立 土壤 沉陷 物理 模型 ,利用 沉陷 剖面 不 同 应 力 
区 土壤 容重 数据 建立 HYDRUS-2D 水 盐 模 型 ,研究 土壤 蒸发 过 程 中 ,土壤 水 分 .总 盐分 、 不 同 溶质 离 
子 在 不 同 沉陷 应 力 区 随 土壤 深度 和 运 移 时 间 的 变化 规律 。 结 果 表 明 :(1) 0~40 cm 深度 ,沉陷 拉 张 


生 用 显著 增强 土壤 燕 发 作用 ,造成 挤 压 区 土壤 含水 率 显 著 大 于 拉 张 区 ,土壤 含水 率 随时 间 变 化 曲 
线 呈 现 典型 的 蒸发 过 程 “ 三 阶段 ”模式 (2) 沉陷 组 各 应 力 区 土壤 总 含 盐 量 不 但 随 深 度 呈 现 强烈 变 
异 ,具体 表现 为 20~40 cm ,60-80 cm 土壤 总 含 盐 量 相对 积聚 ;而 且 右 侧 拉 张 区 总 含 盐 量 积 聚 深度 有 
向 下 迁移 的 趋势 (3) 沉陷 拉 张 区 Ca” SOF Mg" CT 离子 积聚 浓度 大 于 挤 压 区 。Ca SO, CO 士 
壤 离 子 浓 度 随 深度 呈 单 峰 积聚 ,沉陷 拉 张 区 积聚 深度 均 显 著 下 移 。(4) HYDRUS-2D 水 盐 模 型 可 以 


较 好 地 模拟 开采 沉陷 土壤 剖面 中 的 土壤 水 盐 运 移 。 土 壤 含 盐 量 、 含 水 率 的 实测 值 与 模拟 值 率 定 结 
果 为 平均 相对 误差 (ME)<0.5 , 均 方 根 误差 (RMSE )<0.5 ,决定 系数 ( 尼 )>0.95 ,模拟 精度 满足 试验 要 
求 。 研 究 结果 可 为 科学 开展 开采 沉陷 区 生态 修复 工程 、 提 高 生态 自修 复 能 力 提供 重要 理论 依据 。 
X 键 词 : 开采 沉陷 ; 土壤 水 分 运 移 ; 土壤 溶质 运 移 ; HYDRUS-2D 


文章 编号 : 1000 - 6060(2023)09 - 1481 - 12(1481 ~ 1492) 


煤 灵 资源 长 期 占据 我 国 能 源 主体 地 位 。 中 国 
山西 -陕西 -内 蒙古 交界 地 长 期 高 强度 煤 误 开 采 , 造 
成 地 表 沉 陷 ,发 育 不 同 规模 的 腹 颖 ,土壤 内 部 发 育 
微 裂 险 和 陷落 结构 。 由 于 土壤 结构 破坏 ,导致 地 表 
蒸发 量 增 大 ,地 表 水 也 会 沿 裂 颖 癌 深 部 渗 漏 ,造成 
地 表 水 资源 流失 。 同 时 由 于 “ 盐 随 水 走 ” ,地表 沉陷 
不 仅 造 成 土壤 水 分 空间 变异 ,而 且 影 响 土壤 盐分 运 
移 , 加 大 了 土壤 盐 碱 化 风险 。 开 及 沉陷 对 土壤 水 盐 
运 移 的 改变 会 对 土壤 微生物 群落 代谢 活动 地表 植 
被 和 土壤 养分 循环 等 地 表 生 态 环 境 产 生 巨 大 影响 ， 
是 山西 -陕西 -内 蒙古 交界 地 矿区 地 表 生 态 退 化 的 
重要 因素 。 因 此 ,研究 开采 沉陷 地 土壤 水 盐 运 移 
规律 是 矿区 生态 修复 的 重要 组 成 部 分 ,对 科学 评价 
开采 沉陷 的 环境 效应 具有 重要 意义 。 


收 稿 日 期 : 2022-11-20; 修订 日 期 : 2023-01-15 


近年 来 ,国内 外 学 者 “利用 田间 试验 ,物理 模 
型 .数值 模拟 等 方法 研究 采 煤 沉陷 地 土壤 水 盐 运 移 
规律 ,主要 集中 于 沉陷 裂 颖 对 土壤 水 盐 和 植被 根系 
生长 的 影响 .沉陷 盆地 土壤 水 盐 时 空 变异 特征 等 方 
面 。 其 中 ,Tripathi 等 所 以 印度 热带 干旱 地 区 煤田 为 
研究 对 象 ,通过 对 沉陷 斜坡 沉陷 洼地 、 相 邻 未 受 损 
地 区 的 对 比 研 究 发 现 ,沉陷 导致 土壤 容重 土壤 持 
水 能 力 在 沉陷 洼地 增加 ,而 在 沉陷 斜坡 降低 。Xiao 
等 二 以 山东 省 济宁 市 东 滩 煤矿 为 研究 对 象 ,对 土壤 
含水 量 的 空间 变异 性 和 相关 性 进行 分 析 评 价 ,表明 
土壤 含水 量 和 容重 具有 强烈 的 空间 依赖 性 。 张 健 “ 
通过 室内 物理 试验 与 数值 模拟 结合 的 方式 模拟 不 
同 尺度 地 裂缝 周围 土壤 水 分 空间 分 布 变 化 ,研究 结 
果 表 明 地 裂 颖 影响 周围 土壤 水 分 空间 分 布 , 加 速 裂 
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颖 周围 土壤 水 分 散失 , 且 地 裂 颖 尺度 越 大 ,其 影响 
范围 越 大 ; 越 靠近 裂 颖 壁 ,土壤 含水 量 越 低 ; 随 土壤 
深度 增加 ,裂缝 的 影响 作用 逐渐 降低 。 李 走 走 ” 选 
择 典 型 黄土 矿区 大 佛寺 煤矿 为 研究 区 ,采用 遥感 影 
像 反 演 .实地 监测 .光谱 数据 建 模 等 方法 ,从 不 同 尺 
度 分 析 采 煤 沉陷 区 土壤 湿度 的 时 空 演变 特征 ,研究 
结果 表明 采 煤 沉陷 引起 地 表 拉 伸 变 形 和 裂 颖 发 育 ， 
造成 土壤 密实 度 降 低 .土壤 水 分 人 渗 和 蒸发 作用 
增强 。 

随 着 土壤 水 盐 运 移 理 论 模 型 的 深入 研究 以 及 
大 量 田 间 试 验 数据 的 文 撑 ,前 人 开展 了 大 量 水 盐 运 
动 过 程 定量 化 研究 。HYDRUS 水 盐 运 移 模 拟 模型 
是 目前 应 用 最 广泛 的 定量 水 盐 运 移 方法 , 因 其 边界 
条 件 灵活 ,主要 适用 于 室内 土 柱 模拟 8*” 。 其 中 ,Liu 
等 "0 为 了 提高 暗 管 排水 中 HYDRUS-2D 模 拟 的 准确 
性 ,提出 阻力 层 、 实 际 入 渗 面 和 等 效 半径 处 理 3 种 方 
法 ,来 代替 HYDRUS-2D 地 下 管道 的 标准 渗流 边界 ， 
大 大 提高 了 土壤 含 盐 量 含水 率 的 模拟 精度 。 王 世 
明 等 中 通过 HYDRUS 模 拟 防 护林 建设 初期 不 同 矿 
化 度 地 下 水 灌溉 条 件 下 土壤 水 盐分 布 特征 和 时 间 
变化 ,结果 表明 表层 (0~30 cm) 土 壤 含水 量 和 含 盐 
量 受 灌溉 影响 大 ,数值 波动 剧烈 , 且 土 壤 盐分 表 聚 
强烈 ;深层 (50~150 cm) 则 受 洪 溉 的 影响 小 ,数值 波 
动 小 。 于 洋 等 "通过 CMS 和 HYDRUS 的 结合 ,构建 
包 气 带 - 饱 和 带 的 水 盐 耦 合 模拟 模型 ,以 半 王 旱地 
区 宁夏 中 卫 某 速生 林 高 阶地 为 例 ,采用 历史 数据 进 
行 了 验证 ,结果 表明 该 模型 可 有 效 模 拟 包 气 带 - 饱 
和 带 的 水 流 运动 和 和 毛 离子 的 运 移 情 况 。 邹 慧 于 结 
合 室内 试验 HYDRUS-2D 水 盐 模 型 对 神 东 矿 区 煤炭 
开采 中 土壤 水 分 动态 .土壤 水 分 人 渗 .蒸发 及 其 对 
土壤 物理 性 质 的 影响 进行 研究 ,结果 表明 受 地 形 、 
植被 和 开采 的 影响 ,表层 土壤 水 分 呈现 差异 显著 的 
斑 块 状 分 布 格局 ,沉陷 和 裂缝 引起 了 土壤 含水 量 等 
土壤 物理 性 质 的 变化 。HYDRUS-2D 可 以 模拟 土壤 
水 盐 在 水 平方 向 和 任意 复杂 地 形 条 件 下 的 迁移 问 
题 。 因 此 ,对 于 有 效 预 测 点 源 ` 面 源 程度 上 对 陕西 - 
山西 -内 蒙古 交界 地 的 采 煤 沉陷 地 土壤 水 盐 运 移 研 
究 具 有 重要 意义 。 

在 前 人 研究 的 基础 上 ,为 了 揭示 开采 沉陷 微 裂 
际 对 土壤 水 分 和 溶质 运 移 的 影响 ,本 文采 用 不 同 沉 
陷 应 力 区 土壤 容重 的 变化 来 描述 沉陷 微 裂 际 的 变 
化 ,利用 自主 研发 的 土壤 沉陷 模型 进行 模拟 试验 ， 


研究 蒸发 过 程 中 土壤 水 分 .总 盐分 .不 同 洲 质 离子 
在 不 同 沉陷 应 力 区 随 土 壤 深度 和 时 间 变 化 的 规律 ; 
并 通过 采集 沉陷 剖面 不 同 应 力 区 土壤 容重 数据 , 建 
立 基 于 HYDRUS-2D 的 不 同 沉 陷 应 力 区 土壤 水 盐 运 
移 数值 模型 ,为 系统 揭示 开采 沉陷 对 土壤 水 盐 运 移 
的 影响 机 制 、 科 学 评价 开采 沉陷 的 环境 效应 提供 理 
论 依据 。 


1 材料 与 方法 


1.1 沉陷 物理 模型 建立 

本 人 研究 利用 课题 组 自行 设计 的 沉陷 模拟 装置 
进行 沉陷 模拟 试验 (图 1a)。 该 装置 为 一 沿 口 长 方 
体 有 机 玻璃 ,尺寸 :150 cmx30 cmx120 cm。A 底 面 
由 15 个 活动 面板 组 成 。 每 个 活动 面板 尺寸 :10 cmx 
30 emx1 cm。 每 个 活动 面板 由 B 底 面 4 个 活动 螺钉 
文 撑 ,可 通过 扭 动 螺钉 来 调节 活动 面板 升降 高 度 ， 
最 大 可 调 高 度 15 em。 在 A 底面 活动 面板 上 方 分 层 
填充 土壤 ,通过 调 市 螺丝 高 度 模拟 沉陷 过 程 ,A 底面 
两 端 各 留 3 个 活动 面板 ,高 度 不 变 (高 度 为 15 cm)。 
从 左边 第 4 个 活动 面板 回 右 依次 调节 9 个 活动 面板 
高 度 ,模拟 治 谢 面 走向 沉陷 盆地 依次 向 右 推进 ,上 
履 土 层 从 左 至 右 依次 沉陷 JERAM. ERÄ 
面 中 不 同 沉陷 应 力 区 的 划分 是 基于 开采 沉陷 盆地 
的 理论 模型 (图 1b) 而 来 。 沉 陷 盆 地 中 心 部 位 分 布 
挤 压 应 力 ,定义 为 挤 压 区 。 该 挤 压 区 两 侧 分 布 拉 张 
应 力 ,开采 裂缝 集中 发 育 于 该 部 位 ,定义 为 拉 张 
区 。 和 靠近 模型 边界 部 位 受到 物理 模型 边界 效应 影 
响 较 大 ,因此 该 拉 张 区 域 数 据 未 采用 ,定义 为 边 
缘 区 。 
1.2 试验 过 程 和 数据 采集 

试验 在 内 蒙古 科技 大 学 微生物 复 晨 试验 室 进 
行 。 试 验 设 2 种 处 理 , 试 验 组 为 沉陷 组 ,对 照 组 为 无 
沉陷 处 理 , 每 种 处 理 设 3 个 重复 。 试 验 所 用 土壤 是 
采 自 内 蒙古 包头 地 区 的 砂 土 , 过 2 mm 和 去 除 杂 物 
后 ,风干 、 备 用。 根据 粒度 分 析 数 据 , 该 砂 土 与 山 
西 -陕西 -内 蒙古 交界 地 采 煤 沉陷 地 的 风 积 砂粒 度 
组 成 相近 ( 表 1)。 

模型 搭建 :利用 图 1 所 述 沉陷 模型 ,为 避免 无 文 
撑 条 件 下 降低 活动 面板 造成 土 层 垮塌 的 问题 ,首先 
在 A 底面 活动 面板 上 铺设 0.95 cm 厚度 的 石膏 板 ,为 
上 覆土 体 提供 近似 于 覆 宕 的 支撑 力 。 然 后 将 供 试 
土壤 与 矿 化 度 3 gL "(离子 浓度 比 : Na':Cl:CO3 :Ca 
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(a) 不 同 沉陷 应 力 区 特征 模拟 装置 


水 平 位 置 /cm 
0 30 60 90 120 150 
土 层 
位 置 
B 
其 
K 
e 
二 
A 底 面 
螺丝 
t 20 cm 
沉陷 盆地 推进 方向 B 底 


观测 点 


(b) 开采 沉陷 盆地 理论 示意 图 


SS 一 / 


RSK 开采 推进 方向 


底面 为 可 升降 活动 面板 ,A 底 面 上 方 铺设 石膏 板 ;zw os ws ws 分 别 为 随 开 采掘 进 上 有 禾 地 层 沉 陷 盆 地 发 育 过 程 ;xwx 为 最 终 静 态 盆地 。 
图 1 不 同 沉 陷 应 力 区 特征 模拟 装置 及 开采 沉陷 盆地 理论 示意 图 


Fig. 1 Simulation device of different subsidence stress regions characteristics and theoretical diagram of mining subsidence basin 


表 1 供 试 土壤 基本 理化 性 质 
Tab.1 Basic physical and chemical properties of the tested soil 


No 粒 径 占 比 /% 

测试 项 目 一 -一 
粗 砂 (497~2000 jun) 中 砂 (240~497 um) 细 砂 (72~240 jum) 黏土 (<72 um) 

矿区 土壤 16.35 42.03 17.64 

供 试 土壤 16.50 43.34 11.10 


SO7:Mg”=8:4:1:1:1:1) 的 水 溶液 充分 混合 ,维持 最 大 
饱和 含水 率 为 21%。 分 层 装 土 ,每 层 厚 10 cm, 控制 
试验 土壤 容重 为 1.6 gcm^ ,层面 打 毛 ,防止 分 层 效 
应 。 同 时 安装 25 个 监测 点 的 土壤 水 分 采集 器 。 静 
置 24 nh, 待 水 分 平衡 后 ,对 试验 组 进行 沉陷 模拟 。 

沉陷 模拟 过 程 :从 左边 第 4 个 活动 面板 开始 向 
右 依次 下 降 9 个 活动 面板 ,为 达到 最 佳 沉陷 效果 ,每 
隔 1h 从 左 至 右 水 平 推进 10 cm, 即 下 降 1 个 活动 面 
板 。24h 后 达到 稳 沉 阶段 (图 2)。 

蒸发 过 程 实时 监测 :利用 预先 安装 的 EM50 数 
JR de 5TE 传感器 进行 土壤 水 分 和 电导 率 数 据 
实时 监测 。EM50 监测 仪 ( 壤 博 士 RS485 ) 测 量 间隔 
时 间 为 30 min, 仪 妖精 度 : 电 导 率 0~10000 uS m! 
3% FS.m 含水 率 0~50%+2%。 受 5TE 传感器 监测 
范围 限制 ,土壤 含水 率 达 到 12% , 试验 结束 。 蒸 发 
过 程 为 自然 蒸发 ,室内 温度 20~25 % ,湿度 13% ~ 
18%。 

考虑 物理 模型 的 边界 影响 ,本 试验 2 个 边缘 区 


wn NIU 


2 稳 沉 阶 段 模拟 土 柱 实物 
Fig. 2 Physical diagram of simulated soil column 


in steady settlement stage 


数据 未 采用 。 

由 于 传感器 所 测 电导 率 受 土壤 含水 率 影响 ， 
此 采用 土 水 比 1:5 的 土壤 浸 提 液 分 别 测定 25 个 观测 
点 的 土壤 浸 提 液 电 导 率 (ECs,dS.m-)。 土 壤 自动 
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传感器 所 测 电 导 率 (ECsw,dS.m'!) 与 舍 水 率 (0,%) 
的 转换 关系 为 


EC, 70.1120 - 0.12 (1) 
土壤 EC 转换 为 ECis 的 换算 关系 为 : 
ECswEC,s= 22.0190 + 2.82 (2) 


土壤 EC 转换 为 土壤 含 盐 量 (SSC,g'kg) 的 计 
算 公 式 为 “. 
SSC=3.77ECis-0.24 (3) 
13 HYDRUS-2D 水 盐 模型 建立 
1.3.1 土壤 水 分 运动 基本 方程 在 采 煤 沉陷 条 件 
下 ,土壤 水 分 运动 简化 成 径 向 和 垂直 2 个 方向 ,属于 
二 维 运动 水 分 和 人 渗 问 题 。 假 设 试验 土壤 均 质 ,各 问 
同性 。 在 考虑 非 人 饱和 土壤 水 滞后 效应 情况 下 ,土壤 
水 分 运动 规律 的 数学 模型 “” 为: 


00. à ôP., ð op | [9K(O)] 
2r as Ko Bt) d T (4) 


式 中 : 0 为 土壤 体积 含水 率 (ecm em); pH ERK 
负 压 (cm) ; K(9) Jj 3E TRURIL SE S 7K S (em d.) ;2 H 
EAER, EDIE sx DU AER st DIET TR] 

1.3.2 土壤 溶质 运 移 模型 在 均匀 介质 中 ,用 可 控 
的 对 流 弥 散 方程 模拟 非 反应 离子 运 移 ,其 数学 模 
AY USJ; 


age ð (gs xe) -SE 
ot 2 (o D, OX OX (5) 


式 中 ;iy 为 xz 轴 坐标 ;c 为 液体 中 的 溶质 浓度 (g:cm); 
4 为 土 体 中 水 的 流速 ; D; 为 水 动力 弥散 系数 ;i 为 
时 间 。 

1.3.3 土壤 水 分 特征 方程 土壤 的 水 力 特征 方程 一 
般 采 用 VG 模型 ,其 表达 式 " 为 : 


(a) 试 组 验 


100 cm 


自由 排水 边界 


_0h)-0. _ m 
6.7, ag ^( * N (6) 
K(0) "ka -(1- 2 (7) 


式 中 : K, 为 土壤 饱和 导 水 率 (cm*d); 9. 为 土壤 相对 饱 
ME lemem); 0, 为 土壤 相对 饱和 导 水 率 (cmom ) ; 
0. 为 土壤 饱和 体积 含水 率 (cm * em) ;为 负 压 水 头 
(em) :Qn 为 试验 测定 经 验 函 数 ;m=1-1/n;l 为 形状 
系数 ,通常 取 平 均值 0.5。 
1.4 模型 定 解 条 件 

土壤 水 分 运 移 的 边界 条 件 如 图 3 所 示 。 

土壤 溶质 运 移 的 初始 条 件 : 


C(z,0)= C,(z OxzxZ (8) 
土壤 盐分 运动 上 边界 为 大 气 边界 : 
E uu +9C=gC0 (9) 
下 边界 为 浓度 边界 : 
C(zD=CO (10) 


式 中 : C HERRE Elg ke); JEER, 
向 上 为 正 ; q, HEARKE; C, 为 上 边界 电导 率 
IKE C(g-em?); C, 为 下 边界 电导 率 浓 度 (g.cm3); 
为 时 间 。 

15 模型 参数 

1.5.1 土壤 水 力 特 征 参 数 采用 土壤 容重 指标 量化 
地 表 变 形 程度 。 利 用 供 试 土壤 粒 径 ,沉陷 剖面 不 同 
应 力 区 25 处 观测 点 实测 容重 数据 ,由 HYDRUS-2D 
的 Rosetta 程序 自动 生成 土壤 水 力 特 征 参 数 ( 表 2)。 
1.5.2 单元 划分 及 离散 根据 试验 过 程 中 土 柱 沉降 
变化 与 实测 土壤 容重 数据 建立 二 维 模型 ,模拟 沉陷 


(b) 对 照 组 


100 cm 


自由 排水 边界 


图 3 不 同 沉陷 应 力 区 土壤 水 分 和 溶质 运 移 数 学 模型 


Fig.3 Mathematical models of soil moisture and solute transport in different subsidence stress regions 
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表 2 土壤 水 力 特征 参数 


Tab.2 Soil hydraulic characteristic parameters 


处 理 组 0./cm? cm 0./cm - cm? 
左 侧 边缘 区 0.0436 0.3956 
拉 张 B 区 0.0415 0.3607 
挤 压 C 区 0.0413 0.3579 
拉 张 D 区 0.0368 0.2802 
右 侧 边 缘 区 0.0405 0.3468 
对 照 组 0.0375 0.3364 


a n 天 /cm'd 
0.0363 1.6585 96.68 
0.0395 1.6303 61.54 
0.0396 1.6243 59.24 
0.0498 1.4600 17.38 
0.0414 1.5969 50.64 
0.0469 1.6164 53.80 


注 :0. 为 土壤 相对 饱和 导 水 率 ;0. 为 土壤 饱和 体积 含水 率 ;an 为 试验 测定 经 验 函数 ;K. 为 土壤 饱和 导 水 率 。 


土壤 剖面 。 试 验 组 与 对 照 组 模拟 区 域 与 尺寸 如 图 3 
所 示 。 采 用 不 等 间距 网 格 对 模拟 区 域 进行 离散 
化 。 试 验 组 生成 节点 1147 个 ,离散 网 格 2500 个 。 
对 照 组 生成 节点 196 个 ,离散 网 格 4802 个 。 时 间 离 
散 单 位 为 d, 模 拟 时 段 为 2020 年 9 月 18 日 一 2021 年 
11 月 4 日 ,共计 78 d。 设 定 初始 时 间 步 长 0.1 d, 最 小 
时 间 步 长 0.001 do 


2 结果 与 分 析 


2.1 HYDRUS-2D 水 盐 模型 精度 检验 
为 了 确保 HYDRUS-2D 水 盐 模型 的 可 靠 性 ,将 
LE dde 含水 率 实 测 值 与 模型 的 模拟 值 进行 平 
均 相 对 误差 (ME) 、 均 方 根 误差 (RMSE) 和 决定 系数 
(也 ) 检 验 ( 表 3)。 当 ME<0.5、RMSE<0.5、R2>0.95 
时 ,模拟 精度 满足 试验 要 求 。 就 整体 而 言 ,土壤 水 
盐 模拟 值 与 实测 值 吻 合 度 较 好 ,模拟 结果 能 较 好 地 
反映 土壤 水 盐 的 动态 变化 。 这 证 明 HYDRUS-2D 水 


盐 模型 可 以 较 好 地 模拟 不 同 沉 陷 应 力 区 土壤 剖面 
中 的 水 盐 运 移 。 
2.2 不 同 沉陷 应 力 区 土壤 含水 率 时 空 变化 特征 
沉陷 拉 张 作用 增强 了 土壤 蒸发 作用 , 拉 张 区 土 
壤 含 水 率 显著 低 于 挤 压 区 ,土壤 含水 率 随 深度 变化 
显著 (图 4)。<40 cm 深度 土壤 含水 率 随 时 间 变 化 曲 
线 呈 现 典型 燕 发 过 程 “ 三 阶段 " :稳定 蒸发 阶段 (0~5 d), 
蒸发 速率 下 降 阶 段 (5~40 d) RIZE AKA HI BET EE 
(40 4 以 后 )。 在 葵 发 全 过 程 中 ,0~20 cm 沉陷 挤 压 C 
区 土壤 含水 率 显著 高 于 两 侧 沉 陷 拉 张 区 ,与 对 照 区 
差异 较 小 。 这 说 明 沉 陷 前 面 2 个 拉 张 区 土壤 微 裂 除 
发 育 , 增 大 了 土壤 与 大 气 的 接触 面 , 薪 发 作用 强 
烈 。 但 是 ,20~40 cm 、40~60 em 左 侧 拉 张 B 区 和 挤 压 
C 区 的 土壤 含水 率 无 显著 差异 。 类 比 随 开采 推进 方 回 
沉陷 盆地 的 发 育 过 程 ,这 是 由 于 随 着 沉陷 向 右 推进 ， 
左 侧 20~60 em 沉陷 拉 张 裂 颖 闭合 的 原因 。>40 cm 
深度 土壤 含水 率 在 蒸发 作用 初期 显著 增 大 ,可 能 是 
上 层 土 壤 水 分 癌 下 运 移 造成 的 下 层 土 壤 含 水 率 波 


表 3 模型 精度 检验 


Tab.3 Tests of model accuracy 


处 理 组 土 层 深度 /em 人 
RMSE ME R RMSE ME R 

试验 组 0-20 0.03 0.01 0.99 0.04 —0.01 0.95 
20-40 0.03 0.01 0.99 0.04 —0.04 0.96 
40-60 0.02 —0.01 0.98 0.03 0.01 0.95 
60-80 0.05 0.05 0.93 0.03 0.01 0.94 
80-100 0.02 0.04 0.94 0.04 —0.01 0.96 

对 照 组 0-20 0.02 -0.02 0.95 0.03 0.04 0.96 
20-40 0.03 —0.02 0.94 0.04 0.03 0.97 
40-60 0.01 0.01 0.92 0.01 0.03 0.96 
60-80 0.02 —0.02 0.96 0.01 0.03 0.95 
80-100 0.03 —0.03 0.97 0.01 0.02 0.97 


IE: ME 为 平均 相对 误差 ;RMSE 为 均 方 根 误差 ;R 为 决定 系数 。 
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m 拉 张 B 区 
一 一 拉 张 B 区 


e HECK 
---: 挤 压 C 区 


(a) 0~20 cm 


土壤 含水 率 /% 


土壤 含水 率 /% 


时 间 /d 
(e) 80-100 cm 
23 
22 
21 
20 
19 
18 
17 
16 
15 


土壤 含水 率 /% 


A 拉 张 D 区 
iiis 拉 张 D 区 


(b) 20~40 cm 
2r I; II 


v 对 照 组 
对 照 组 


士 壤 含水 率 /% 
3 


土壤 含水 率 /% 


时 间 /d 


时 间 /d 


注 : 工区 域 为 稳定 蒸发 阶段 (0~5 d); 工区 域 为 蒸发 速率 下 降 阶 段 ($~40 d); 亚 区 域 为 莹 发 水 分 扩散 阶段 (40 d 以 后 )。 
图 4 不 同 深度 不 同 沉陷 应 力 区 土壤 含水 率 随 时 间 变 化 


Fig. 4 Variation of soil moisture with time in different subsidence stress regions at different depths 


动 。>60 cm 深度 挤 压 C 区 土壤 含水 率 仍 然 显著 大 
于 两 侧 拉 张 区 。 这 说 明 沉陷 剖面 深部 拉 张 区 对 土 
壤 含 水 率 的 扰动 时 效 性 比 浅 层 更 久 。 

蒸发 试验 结束 时 ,不 同 沉陷 应 力 区 土壤 含水 率 
随 深 度 的 分 布 剖 面 图 (图 5) 表 明 ,0~40 em 深度 , 沉 
陷 拉 张 作用 显著 增强 土壤 蒸发 作用 ,造成 挤 压 C 区 
土壤 含水 率 显 著 大 于 拉 张 区 。40~60 cm 深度 不 同 
沉陷 应 力 区 土壤 含水 率 无 显著 差异 。>60 cm 深度 
土壤 含水 率 又 呈现 随 沉陷 应 力 区 不 同 而 分 异 的 
特征 。 


2.3 不 同 沉 陷 应 力 区 土壤 总 含 盐 量 时 空 变化 特征 

不 同 深度 土壤 总 含 盐 量 随时 间 逐 渐 减 小 ,以 0~ 
20 cm 深度 土壤 随时 间 脱 盐 量 最 大 (图 6)。 这 是 表 
层 土壤 强烈 蒸发 作用 导致 盐分 在 0~5 cm 深度 形成 
盐 壳 的 结果 。0~60 cm 深度 左 侧 沉 陷 拉 张 B 区 
土壤 脱盐 量 随时 间 逐 渐 增 大 ,与 挤 压 C 区 差异 显 
著 。>60 cm 深度 土壤 总 含 盐 量 随时 间 变 化 曲线 无 
显著 差异 。 

沉陷 不 同 应 力 区 土壤 总 含 盐 量 随 深度 
变异 (图 7)。 蒸 发 试验 10 d,0~20 cm 土壤 盐分 
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(a) 试验 组 (b) 对 照 组 
水 平 位 置 /cm 水 平 位 置 /cm 
0 30 60 90 120 150 0 30... 60 . 90 


土壤 含水 量 /cm? - cm? THES Kien? : cm? 


土壤 深度 /cm 
g 


80 
100 


图 5 第 78 d 时 土壤 水 分 剖面 模拟 


Fig. 5 Simulation of soil moisture profile on the 78" day 


m 拉 张 B 区 © PECK 4 拉 张 D 区 v 对 照 组 
一 一 拉 张 B 区 ---: 挤 压 C 区 ..…. 拉 张 D 区 对 照 组 


(a) 0-20 cm (b) 20-40 cm 


土壤 总 含 盐 量 /gkg! 
土壤 总 舍 盐 量 /g*kg! 


土壤 总 舍 盐 量 /gkg 1 
THEE E Eb e kg! 


时 间 /d 时 间 /d 
(e) 80~100 cm 


E 
i 
E 
g 
n] 
EK 
-H 
0 20 40 60 80 
时 间 /d 


图 6 不 同 深度 不 同 沉陷 应 力 区 土壤 总 含 盐 量 变化 


Fig. 6 Changes of soil total salt content in different subsidence stress regions at different depths 
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(a) 10 d 试 验 组 水 平 位 置 /cm (b) 10 d 对 照 组 水 平 位 置 /cm 
30 60 90 120 — 150 30 60 90 120 — 150 
0 土壤 总 含 盐 量 /gker 0 HER Ard lg kg 
. 4.3 
20 20 — 3.9 
g e E —EH 
9 40 : 9 40 pi 
EN : EA — HX 
S oo S 60 m, 
Ei : Ek um ^ 
H 1.8 H FE; 
80 1.4 80 = 1.4 
1.0 — $E 
0.7 0.7 
100 03 100 Em, 
(c) 50 d 试 验 组 水 平 位 置 /cm (d) 50 d 对 照 组 水 平 位 置 /cm 
30 60 90 120 — 150 30 60 90 120 — 150 
0 土壤 总 含 盐 量 /gkg" 0 土壤 总 合 盐 量 /g*kg” 
4.3 4.3 
20 m., 20 m. 
m". LE 
g -一 g — 
En ——pX E E 
. 2.5 
xe " [a 21 K 60 = 2.1 
gl m. S m. 
H —m H E. 
一 二 —Ix 
100 03 100 03 
(e) 78 d 试 验 组 水 平 位 置 /cm (f) 78 d 对 照 组 水 平 位 置 /cm 
30 60 90 120 — 150 30 60 90 120 — 150 
0 HAE RU kg 0 土壤 总 含 盐 量 /g"kg-! 
4.3 4.3 
20 3.9 20 E 
8 3.6 B —IY. 
2 ə — E 
其 40 F i 40 一 了 
EE 25 ES —— px: 
i 60 24 ii 60 E. 
r - 
80 LO 80 E, 
0.7 — 
100 03 100 mE. 


图 7 不 同 沉陷 应 力 区 土壤 总 含 盐 量 模拟 剖面 


Fig. 7. Simulation profile of soil total salt content in different subsidence stress regions 


发 作用 影响 , 回 土 壤 最 表层 迁移 ,沉陷 土壤 总 售 盐 
量 高 于 对 照 组 , 右 侧 拉 张 D 区 盐分 积累 量 显 著 高 于 
其 他 应 力 区 ,盐分 最 终 在 土壤 表层 0~5 cm 形成 盐 
这。 在 50 d、78 d, 对 比 对 照 组 ,沉陷 组 各 应 力 区 土 
壤 总 含 盐 量 不 但 随 深度 呈现 强烈 变异 ,具体 表现 为 
20-40 cm 、60~80 cm 土壤 总 含 盐 量 相对 积聚 ;而且 
右 侧 拉 张 区 总 含 盐 量 积 聚 次 度 有 向 下 迁移 的 趋 
势 。 这 表明 沉陷 区 深部 土壤 水 分 不 均匀 分 布 会 导 
致 土壤 盐分 随 深度 的 变异 , 拉 张 区 土壤 盐分 向 下 迁 
移 ,减缓 了 沉陷 地 盐 碱 化 风险 。 
2.4 不 同 沉 陷 应 力 区 6 种 土壤 溶质 离子 分 布 特征 
沉陷 拉 张 区 和 挤 压 区 显著 改变 了 土壤 溶质 离 


的 积聚 效应 (图 8)。 在 对 照 组 中 ,Ca*、SO7、CO3 土 
壤 离 子 浓 度 均 随 深 度 呈 现 典型 双 峰 式 积聚 效应 , 积 
聚 深 度 为 20~40 cm .60~80 cm。 对 比 对 照 组 ,沉陷 
MIK Ca” SOF, COF 土壤 离子 浓度 随 次 度 呈 单 峰 
积聚 , 且 Ca SO 离子 在 沉陷 拉 张 区 积聚 深度 下 移 
至 40~60 cm, 挤 压 C 区 无 显著 积聚 效应 ;CO? 离子 在 
沉陷 挤 压 C 区 和 拉 张 区 具有 相似 分 布 规律 , 仅 在 
60-80 cm 积聚 ,呈现 显著 向 下 迁移 的 趋势 。 在 对 照 
组 中 ,Me ,Na',CE 土壤 离子 浓度 以 单 峰 式 积聚 。 对 
比 对 照 组 ,在 沉陷 拉 张 区 3 种 离子 积聚 深度 均 上 移 ; 
在 沉陷 挤 压 区 Mg" , CE 离子 随 深 度 无 显著 积聚 效 
应 ;Na 离子 在 沉陷 挤 压 区 呈 双 峰 式 积聚 。 沉 陷 拉 张 


子 随 次 度 的 变化 ,不同 的 涂 质 离子 随 深 度 呈 现 不 同 


[X Ca” SOF Mg" CT 离子 积聚 浓度 大 于 挤 压 区 。 
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8 不 同 沉陷 应 


Fig. 8 Distribution profiles of six solute ions under different subsidence stress 


3 讨论 


3.1 沉陷 应 力 区 分 布 对 土壤 水 分 运 移 的 影响 机 制 
在 中 国 山 西 -陕西 -内 蒙古 三 省 交界 地 ,高 强度 
的 井 工 开采 造成 风 积 沙 地 黄土 覆盖 区 大 面积 沉陷 ”。 
在 开采 沉陷 扰动 下 和 稳定 沉陷 盆地 形成 过 程 中 , 土 
体内 部 形成 2 类 裂 颖 ,一 类 为 边缘 裂缝 ,在 开采 沉陷 
过 程 中 不 断 加 深 , 且 不 会 在 稳 沉 后 消失 。 男 一 类 为 
动态 裂缝 , 随 开采 沉陷 推进 ,呈现 拉 张 -闭合 - 拉 张 
的 动态 变化 过 程 ”。 而 本 文 室内 模拟 试验 中 ,通过 
依次 调 广 活动 面板 模拟 土壤 受 开采 沉陷 应 力作 用 
形成 的 拉 张 区 和 挤 压 区 ,最 终 形 成 稳定 沉陷 盆地 。 
土壤 容重 指标 量化 沉陷 裂缝 分 布 和 发 育 情况 ,不 同 
沉陷 应 力 区 导致 土壤 容重 发 生 显 车 变化 ,试验 研究 
发 现 开采 沉陷 对 土壤 水 盐 运 移 的 主要 影响 因素 是 
土 体 裂 缝 的 发 育 情况 ,包括 宏观 裂缝 和 土 体内 部 的 
微 裂 际 。 这 些微 错 际 为 土壤 水 分 运 移 创造 了 优先 
路 径 ,加 速 土 壤 水 盐 运 移 过 程 。 这 也 侧面 印证 了 Hu 
等 ”的 结论 ,土壤 含水 量 和 容重 具有 强烈 的 空间 依 
赖 性 。 这 成 为 干旱 少雨 光照 强度 大 的 西部 地 区 地 


1.5 


20 2.5 d 02 03 0.4 0.5 


土壤 深度 /cm 


影响 下 6 种 溶质 离子 分 布 剖面 


表 生 态 损伤 的 重要 因素 。 

毕 银 丽 等 2 发 现在 沙 地 沉陷 区 ,开采 沉陷 增强 
了 60~80 cm 土壤 水 分 的 变异 性 ,沉陷 区 土壤 蒸发 量 
均 明 显 大 于 未 沉陷 区 。 本 文 室内 模拟 试验 结果 表 
HH ,0~40 cm 深度 ,沉陷 拉 张 作用 显著 增强 土壤 蒸 
发 ,造成 挤 压 区 土壤 含水 率 显 著 大 于 拉 张 区 。 这 说 
明 沉 陷 拉 张 应力 导 致 土壤 内 部 微 裂 际 发 育 , 增 加 土 
壤 与 大 气 的 接触 面积 ,蒸发 作用 增强 ,加 速 土壤 水 
分 运 移 。 本 试验 中 ,>60 em 深度 土壤 含水 率 又 呈现 
随 沉陷 应 力 区 不 同 而 分 异 的 特征 。 综 合 分 析 , 浅 煤 
层 开采 沉陷 会 在 60~80 cm 深度 形成 不 同 程度 垮塌 ， 
垮塌 在 沉陷 土壤 剖面 的 水 平 位 置 具 有 随机 性 ,这 造 
成 土壤 水 分 在 垮塌 区 变异 性 增强 。 

本 试验 研究 结果 还 表明 ,开采 沉陷 对 土壤 水 分 
影响 具有 时 效 性 ,土壤 水 分 在 一 定时 间 后 趋 于 稳定 
状态 。 这 一 规律 表明 开采 沉陷 对 地 表土 壤 、 植 物 根 
系 万 至 地 表 生 态 系 统 的 影响 均 随 时 间 变 化 。 随 着 
动态 裂 颖 的 闭合 ,开采 沉陷 区 地 表 生 态 系统 表现 出 
显著 的 自修 复 特 点 。Druzbicka 等 中 通过 对 干旱 、 半 
干旱 地 区 2 个 废弃 矿山 的 研究 发 现 , 采 矿 引 起 岩层 
暴露 与 沉陷 促进 了 土地 盐 碱 化 ,但 通过 长 时 间 自 然 
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演 替 过 程 , 受 干 扰 土 地 自然 恢复 ,地 区 生物 多 样 性 
增强 。 因 此 ,在 开采 沉陷 地 生态 修复 过 程 中 ,我 们 
必须 充分 认识 沉陷 地 土壤 水 分 空间 变异 性 与 时 效 
性 ,有 机 结合 人 工 修复 技术 ,加快 沉 陷 区 生态 自修 
复 进 程 。 
3.2 开采 沉陷 促进 干旱 地 区 土壤 盐分 向 下 运 移 

“ 盐 随 水 走 ” ,土壤 盐 分 的 运 移 均 随 水 分 的 运动 
而 迁移 后 。 因 此 , 拉 张 应 力 区 土壤 微 裂 辽 的 发 育 不 
仪 增强 土壤 蒸发 作用 ,而 且 势 必 影 响 土 壤 盐 分 运 
移 ,这 是 导致 山西 -陕西 -内 蒙古 三 省 交界 地 土壤 盐 
碱 化 的 重要 因素 。 

Ashish 等 2 发 现 底部 即 沉陷 挤 压 应 力 区 中 心 位 
置 0~60 cm 深度 土壤 电导 率 显 苦 增 加 ,也 证 明 沉 陷 
应 力 增 加 了 土壤 盐分 的 空间 变异 性 。 本 文 试验 研 
究 结 果 表 明 ,沉陷 应 力作 用 对 土壤 含 盐 量 影响 深度 
远大 于 土壤 含水 率 。0~60 cm 左 侧 沉陷 拉 张 区 土壤 
脱盐 量 随 时 间 逐 渐 增 大 ,与 挤 压 区 差异 显著 。 这 说 
明 开 采 沉 陷 对 土壤 盐分 空间 异 质 性 的 影响 比 土壤 
水 分 更 显著 。 史 文 娟 等 ”通过 室内 土 柱 模拟 试验 ， 
研究 蒸发 条 件 下 不 同 层 位 夹 砂 层 土壤 剖面 盐 离 子 
的 动态 变化 情况 ,结果 发 现 盐分 离子 向 上 迁移 的 速 
度 随 时 间 的 增加 而 减 小 。 此 外 ,本 试验 结果 还 表 
HH ,沉陷 应 力 使 0~$ cm 土壤 积 盐 明显 。 在 前 人 研究 
的 基础 上 ,本 文 研究 结果 证 明了 开采 沉陷 作用 下 土壤 
表层 盐 这 的 形成 进一步 阻止 水 分 和 盐分 向 上 迁移 。 

沉陷 应 力作 用 影响 土壤 质量 ,导致 沉陷 区 域内 
土壤 物理 和 化 学 性 质 受 影响 ,尤其 在 中 国 山西 - 陕 
西 -内 蒙古 三 省 交界 地 易 造 成 土地 盐 碱 化 ”。 土 壤 
性 质 受 到 自然 因素 和 人 为 因素 影响 ,具有 极 强 的 空 
间 变 异性 。 本 试验 研究 发 现 ,沉陷 区 深部 土壤 水 分 
不 均匀 分 布 会 导致 土壤 盐分 随 深度 的 变异 , 拉 张 区 
土壤 盐分 向 下 迁移 。 这 与 毕 银 丽 等 ”利用 箱 状 试 
验 模拟 开采 沉陷 裂 颖 ,进行 HYDRUS-2D 数值 模拟 
研究 结果 相似 ,人 研究 结果 表明 裂缝 区 含 盐 量 向 下 迁 
fe ,证 明 开 采 沉 陷 裂 颖 抑制 了 次 生 盐 碱 化 发 生 。 本 
研究 发 现 沉 陷 拉 张 区 Ca” SOF Mge” CL 离 子 积聚 浓 
度 大 于 挤 压 区 ,Ca*、SO? CO 土壤 离子 浓度 随 深度 
呈 单 峰 积 聚 ,沉陷 拉 张 区 积聚 深度 均 显 著 下 移 。 这 
进一步 证 明 不 同 盐分 离子 具有 向 下 运 移 的 趋势 , 且 
离子 类 型 影响 盐分 的 下 移 规律 。 

土壤 养分 的 变异 是 土壤 水 盐 变 化 的 后 效 作用 5 。 
开采 沉陷 应 力作 用 影响 并 改变 土壤 水 盐 运 移 ,最 终 


影响 土壤 养分 ,抑制 植物 生长 。Zhao 等 ”对 黄土 地 
区 采 煤 沉陷 区 采样 分 析 土 壤 有 机 质 与 全 须 ,结果 表 
明 沉 陷 降 低 了 0~20 cm 深度 上 土壤 有 机 质 与 全 损 含 
量 。 地 面 沉降 增 大 了 有 机 质 和 全 毛 的 空间 异 质 
性 。 裂 缝 宽度、 裂缝 深度 和 取样 点 与 裂缝 边缘 的 距 
离 是 导致 土壤 有 机 质 和 全 和 毛 水 平 空 间 异 质 性 的 主 
影响 因素 。 

因此 ,开采 沉陷 对 地 表 生 态 系 统 的 影响 程度 很 
大 。 针 对 开采 沉 隐 不 同 应 力 区 的 生态 恢复 ,要 开展 
“分 区 "生态 修复 模式 ,以 土壤 自修 复 为 基础 ,建立 
差别 化 生态 修复 模式 。 


4 结论 


本 文 聚 焦 不 同 沉陷 应 力 对 沉陷 盆地 水 盐 运 移 
的 影响 ,利用 自主 研发 的 室内 土 柱 模拟 装置 进行 沉 
陷 模 拟 试 验 , 采 用 土壤 容重 指标 量化 沉陷 裂缝 分 布 
和 发 育 情 况 ,研究 蒸发 过 程 中 土壤 水 盐 在 不 同 沉陷 
应 力 区 的 运 移 规律 ,揭示 开采 沉陷 应 力 分 布 对 土壤 
水 盐 运 移 的 影响 机 制 ,得 到 如 下 结论 : 

(1) 开采 沉陷 不 同 应 力作 用 下 ,土壤 含水 率 随 
时 间 变 化 曲线 呈现 典型 蒸发 过 程 “ 三 阶段 ?模式 。 
0-40 em 深度 ,沉陷 拉 张 作用 显著 增强 土壤 蒸发 
作用 。 

(2) 沉陷 组 各 应 力 区 土壤 总 含 盐 量 不 但 随 深度 
呈现 强烈 变异 ,而 且 右 侧 拉 张 区 总 含 盐 量 积聚 深度 
有 向 下 迁移 的 趋势 。 

(3) 沉陷 拉 张 区 Ca、 SOF Mg、CL 离 子 积聚 浓 
EKT ERK., Ca" SOF CO 土壤 离子 浓度 随 深 
度 呈 单 峰 积 聚 ,沉陷 拉 张 区 积聚 深度 均 显 著 下 移 。 

(4) HYDRUS-2D 水 盐 模型 可 以 较 好 模拟 开采 
沉陷 土壤 剖面 中 的 土壤 水 盐 运 移 。 模 拟 精 度 满足 

(5) 在 山西 -陕西 -内 蒙古 交界 地 采 煤 沉陷 区 , 
土壤 盐分 沿 沉陷 拉 张 应 力 区 迁移 的 重要 过 程 可 以 
很 大 程度 上 减缓 因 强 烈 燕 发 作用 造成 的 土地 盐 碱 
化 问题 。 这 对 科学 开展 开采 沉陷 区 生态 修复 工程 、 
建立 差别 化 生态 修复 模式 、 加 速生 态 自修 复 能 力 提 
供 重要 理论 依据 。 
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Abstract: Mining subsidence in the Shanxi- Shaanxi- Inner Mongolia border region of northwestern China has 
caused ecological damage by strongly disturbing the surface soil. To explore the mechanism by which mining sub- 
sidence influences soil water and salt transport, physical soil subsidence and numerical soil-water and salt-trans- 
port models were established in the HYDRUS-2D water-salt model using soil bulk density data from different 
stress areas of the subsidence profile. The variation in soil moisture, total salt, and different solute-ion content 
with soil depth and migration time in different subsidence stress regions during soil evaporation was studied. The 
results showed that: (1) At 0—40 cm depth, the effect of subsidence and tension significantly enhanced soil evapo- 
ration, resulting in soil water content in the extrusion zone significantly higher than the tension zone. (2) The total 
soil salt content in each subsidence stress region varied strongly with depth, and salt accumulation was observed 
in the 20-40 cm and 60-80 cm soil layers. Moreover, the salt accumulation depth in the right tensile area showed 
a downward migration trend. (3) The accumulation concentrations of Ca^, SO? , Mg", and CI in the subsidence 
tensile zone were greater than the extrusion zone. The concentrations of Ca^', SO? , and CO; showed a single 
peak of accumulation with depth and a significant decrease in accumulation depth in the subsidence and exten- 
sion areas. (4) The HYDRUS-2D water-salt model can accurately simulate soil water and salt transport in mining 
subsidence soil profiles. Compared with the measured soil salt and water content values, the average relative er- 
ror, root mean square error, and coefficient of determination of the simulated values were respectively as fol- 
lows: «0.5, «0.5, and 70.95. In the mining subsidence areas of the Shanxi- Shaanxi- Inner Mongolia border re- 
gion, salt migration through the subsidence tensile stress area of the soil can greatly alleviate the problem of land 
salinization caused by strong evaporation. This study thus provides an important theoretical basis for implement- 
ing science-based ecological restoration projects in mining subsidence areas, for establishing differentiated eco- 
logical restoration models, and for accelerating ecological self-healing capabilities. 


Key words: mining subsidence; soil moisture migration; soil solute transport; HYDRUS-2D 


